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1. Introduction 

Fluctuations of the electron spin density play a predominant role in the 
thermodynamics of ferromagnetic metals (see, e.g., [1] and references therein). 
Most of the progress in the spin-fluctuation theory (SFT) has been achieved 
within the functional-integral approach [2, 3]. Using the single-site and static 
approximation, Hubbard [4], Hasegawa [5] and Grebennikov et al. [6] obtained a 
quantitative description of magnetic properties at finite temperatures, which was 
much better than in the Stoner mean- field theory (figure 1, left). To go beyond the 
single-site approximation Hertz and Klenin [7] suggested a self-consistent Gaussian 
approximation. However, it used the static long-wave limit and was restricted to 
paramagnets. To ferromagnetic metals, the variational approach [7] was extended 
by Grebennikov [8] who took into account dynamics and space correlations of the 
fluctuations, the latter again only in the paramagnetic state. A complete dynamic 
non-local approximation to the SFT for ferromagnetic metals was developed by Reser 
and Grebennikov [9]. 

The Gaussian SFT [9] gives a good agreement with experiment over a wide range 
of temperatures. However, at high temperatures the Gaussian approximation yields a 
discontinuous change in magnetic characteristics (see [10] and references therein). The 
main reason for the first-order phase transition is that Gaussian approximation implies 
independent 'harmonic' fluctuations of the spin density and thus fails to account for 
their interaction. 

The first-order phase transition has been observed in various versions of the 
self-consistent renormalization (SCR) theory of spin fluctuations, developed for weak 
ferromagnetic metals (for a review, see [1]). Particularly, in [11] it was argued that the 
first-order discontinuity in the SCR theory can be eliminated by taking into account 
the rotational invariance of the system. This leads to two integro-differential equations 
for the longitudinal X|| and transverse xj- susceptibilities. A simple relation that 
couples x\\ s-iid xi- suggested in [12] from the assumption that the total local spin 
fluctuation, i.e. the sum of the zero-point and thermal spin fluctuations, is conserved, 
as it is in the Heisenberg local moment theory and may be somehow justified for weak 
ferromagnets. 

In the present paper, wo improve the coupling of the fluctuations in the Gaussian 
SFT [9] by taking into account high-order terms of the free energy of electrons F{V) 
in the fluctuating exchange field V (sec [13] for a brief summary). First, in the 
fourth-order Taylor expansion of the free energy F{V), we take a partial average with 
respect to /S.V in the third- and fourth-order terms replacing them by a linear and 
quadratic terms, respectively. Adding these corrections to the Taylor terms of the 
first- and second-order, we come to the extended function F{V). The best quadratic 
approximation is constructed as in [7-9], with the help of the free energy minimum 
principle [14], but using the first- and second-order derivatives of the extended function 
F{V). In the final computational formulae, the third-order term renormalizes the 
mean field, and fourth-order term renormalizes the susceptibility (this includes the 
Gaussian SFT [9] as a special case with the renormalizations set to unity). 

The fundamental difiference between our treatment of the high-order terms and 
the previous ones is that in our approach the ferromagnetic state is changed self- 
consistently (for treatments of the fourth-order term in the paramagnetic state see 
[1,15-17] and references therein). Another advantage of our approximation to the SFT 
is that the integral equation for the mean Green function (coherent potential equation) 
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is finally reduced to a system of four nonlinear equations with four unknowns, which is 
only slightly more complex than the Stoncr mean-field theory. It is significant to note 
that, solving the coherent potential equation directly requires a number of additional 
simplifications, such as neglect of rotational invariance and the mode-mode 'frequency' 
interactions, and, most important, the single-site approximation [18,19]. 

The extended Gaussian approximation of the SFT is applied to numerical 
calculations of magnetic properties of Feo.65Nio.35 Invar alloy at finite temperatures. 
This alloy has been intensively studied recently (see, e.g., [20] and references therein), 
but mostly at zero temperature, i.e. without the quantum statistics. Our choice of the 
Fe-Ni Invar to illustrate the possibilities of the extended SFT is motivated by problems 
of temperature dependence which were found in the quantum-statistical treatment of 
this Invar [10,21-23]. (Obviously, before that the new method has been tested on a 
simplified clean system, such as elemental Fe.) 

It is known that the Fe-Ni Invar is a complex disordered system. However, the 
comparison of the calculation results for the disordered alloy Feo.65Nio.35 [10,21] and 
ordered compound Fe3Ni [22,23] showed that the effect of disorder in the filling of sites 
with Fe and Ni atoms on the magnetic properties of the Fe-Ni Invar is insignificant. 
This conclusion agrees with earlier results for the ordered and disordered Feo.72Pto.28 
Invars (see, e.g., [24], table 10-1). The weak influence of the atomic disorder on 
the magnetic properties of the Fe-Ni Invar at finite temperatures is explained by the 
integral dependence on the electronic energy structure in all the equations of the SFT. 
The details of the initial density of states (DOS) do not exert the decisive effect on 
the results of the calculation. 



2. Quadratic approximation taking into account high-order terms 

2.1. Free energy of electrons in a fluctuating field 

The Stratonovich-Hubbard transformation [2, 3] replaces the pair interaction of 
electrons by the interaction of the electrons with the exchange field§ 

V^{V,,V2,...), Vj = J2v;{T)a", (1) 

a 

fluctuating in space (see figure 1, right) and in 'time' t £ [0, 1/T] (j being the site 
number, cr" the Pauli matrix, a = x,y,z, and T the temperature, in energy units). 
Hence the partition function can be written as a functional integral 

Z = e--^/^ =(^J e-^°(^)/^ DV{t)^ ' J e-^°(^)/^ e"^^'^)/^ DVir), (2) 
where 

is the energy of the exchange field, and 

Fi{V)=TTrlnG{V) (4) 

is the free energy of non-interacting electrons in the field V expressed in terms of the 
Green function 

G{V) = {z + n-Ho-V)-\ (5) 
§ We neglect the charge field, same as in [9]. 
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Figure 1. A schematic picture of exchange field configurations in 
the Stoner mean-field theory and SFT. 

Here A^d = 5 is the number of degenerate d bands, U is the single-site electron-electron 
interaction constant, z is the energy variable, /i is the chemical potential, and Hq is the 
sum of kinetic and potential energy of the non-interacting band electrons. Symbol Sp 
denotes the trace over spin indices, and Tr stands for the full matrix trace, independent 
of the particular representation. To simplify the notation, we omit the band index 
and write the prcfactor in the trace Tr. In expression (4), we omit the terms 
independent of the field V (for details, see Appendix A). 
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2.2. Partial averaging of high order terms of the free energy 

Since Hq and V cannot be diagonalized simultaneously in either coordinatc-'timc' or 
momentum- 'frequency' spaces, exact formulae (4) and (5) are of little use to calculate 
the partition function (2) without an appropriate approximation. Therefore, integral 
(2) is replaced by a Gaussian integral or, equivalently, the exact expression 

FiV) = Fo{V)+Fi{V) = rTr(^^ +lnG(l/)^ (6) 

is replaced by a quadratic form: 

F{V) ^^(2) (Y^^ ^ Tr( AV^^AV^) . AV = V -V. (7) 
Hence, we arrive to the Gaussian fluctuating field V with the probability density 



p{V) oc exp(^-^Tr{AVAAV)^ , 



(8) 



the mean field V and matrix A being the unknown parameters of the SFT to be 
determined. Gompared to conventional Gaussian SFTs, we take into account the 
'anharmonicity' of the fluctuations by renormalizing the parameters V and A. 

First, we expand function (6) in Taylor series to the fourth order in AV = V — V: 

F{V) « F{V) + T Tr + G{V)Av) + It Tr f + {G{V)AV) ' 



UT ' ' J 2 \UT 
+ ^TTv{G{V)AV)'' + ^TTv{G{V)AV)\ (9) 
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where V is some arbitrary value of the exchange field. For the third-order term, we 
define the Gaussian partial averaging by the formula 

Ti{GAVGAVGAV) « Tr{GAVGAyGAV) + Ti{GAVGAVGAy) 

+ TT{GAVGAVGAy) = 3 Tr{GAVGAyGAV) (10) 

and for the fourth-order term by the formula 

Tr (GAVGAVGAVGAV) 

f« ATr{GAVGAyGAVGAV) + 2 Tr {GAVGAVGAVGAV) 

-2 Tr{GAVGAyGAVGAy) - Tr {GA VGAVGAyGAV ) , (11) 

where the underbrace denotes the averaging with the Gaussian probability density (8): 

(. . . A V ... AV ...) = J{...AV...AV.. .)p{V) DV. 

In formulae (10) and (11), the combinatorial multipliers correspond to the number 
of all possible pairings of the terms AV with the cyclic invariance of the trace taken 
into account. In formula (11), the last two terms are chosen so that both sides have 
equal mean values. These terms will be omitted for brevity since they do not contain 
the variable AV and lead only to a change of unimportant free term in expansion (9). 
Furthermore, we approximate the second term in the right-hand side of (11) by the 
first one. Thus, relation (11) takes the form 

Tr {GAVGAVGAVGAV) ^ &TT:{GA yGAV GAVGAV). (12) 

Using formulae (10) and (12), we transform the third- and fourth-order terms in 
expansion (9) into corrections to the first- and second-order terms, respectively: 

fVAV ~ . -. . \ 

F{V) » F{V) -\-TtA + G{V)AV + G{V)A yG{V)AV G{V)AV] (13) 

1 /Ay2 . , , . . . \ 

-F-TTr — — + G{V)AVG{y)AV + 2,G{V)AVG{V)AVG{V)AVG{V)AV . 

^ \. U J. I I / 

2.3. Quadratic approximation based on the free energy minimum principle 

Following [7-9], the best quadratic approximation F^'^^iV) is constructed using the 
free energy minimum principle [14]. 

Define the average with the Gaussian density (8) by the formula 



(...) 



j exp(-F(2)(y)/T) DV^ j{...) exp(-F(2)(y)/T) DV, 
where F^'^^iV) is a quadratic function of the form (7). Then the identity 

exp(-F(y)/r)Dy = / exp(-(F(y)-F(2)(y))/T)exp(-F(2)(y)/T) w 



can be rewritten as 



y^exp(-F(2)(y)/r)£)yj j exp(— F(y)/r) DV 

exp(-(F(y)-F(2)(y))/T)). 
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Applying inequality (exp/) > cxp(/), / being an arbitrary set of real quantities, and 
taking the logarithm, we come to the upper bound for the free energy: 

jr < jr{2) ^ ^p^y^ _ p(2) ^y^-^^ ( 

where 

T = -Tin j exp(-F(y)/T) DV, T'^^) ^ j e^^[-F^'^\V)/T) DV. 

To get the 'best' approximation in the class of all quadratic functions (7), one 
minimizes the right-hand side of (14) to obtain the equation (for details, see [25]) 

so that V is equal to the mean field V = (V^). The matrix of the quadratic form is 
given by 

1 /d'^F{V)\ 



In the present paper, we apply the formulae (15) and (16) to the modified function 
(13) instead of the original free energy (6), as it was done in [9]. In expression (13), 
we average over V with the Gaussian density (8), everywhere but in AV , and replace 
the Ay = V — V terms by AV = V — V . The averaged linear term in (13) annihilates: 

TTv{^^- + {G{V))AV + {G{V)A yG{V)AV G{V))Av) = (17) 

identically over AV . Hence the quadratic form F^^' {V) contains only the second-order 
term: 

F{V) ^ F('\V) = ^TTri^-j^ + {G{V)AVG{V)AV) 

+ 3{G{V)A yG{V)Ay G{V)AVG{V)AV)] . (18) 

To simplify expressions (17) and (18) one step further, we introduce yet another 
partial averaging: 

GA VGAV GAV « (Tr(l))"^ Tv(^GA yGAy ) GAV. 

Finally, wc replace the average over V of the product of the Green functions G{V) by 
the product of mean Green functions {G{V)) = G. Hence relation (17) reduces to 

U-^Tt{VAV) + {1 + tj)TTt{GAV) =0 (19) 

and quadratic form (18) transforms to 

F(2)(y) = J-Tt{AV'^) + ^{1 + 3tj)TTi{GAVGAV), (20) 

where the correction coefficient is 

T] = (Tr(l))"^ Tr(GAyGAyy (21) 
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2.4- Equations for the mean field and chemical potential 

In the ferromagnetic state, we choose the ^;-axis along the direction of the mean field: 
V = V^a\ 

Then the mean Green function is spin-diagonal due to 

= 0, = 0. (22) 

Hence, using the well-known relation for the Pauli matrices: Sp((t'*(t'') = 25a^, we 
rewrite equation (19) as 

[/-I Tr(T/^Ay^) + (l + ?7)TTr(G^Ay^) = 0. (23) 

In order to transform (23) to the mean-field equation, we use the momentum- 
'frequency' representation. Since the exchange field (1) is diagonal in the coordinate- 
'time' representation, its Fourier transform is homogeneous: 

^V^'nn' = ^l?-k', n-n'' (24) 

where k is the wave vector, taking values in the Brillouin zone, and w„ = (2n + V)tiT 
are the thermodynamic 'frequencies'. Furthermore, the mean field is a constant in 
the coordinate-'time' representation, and hence its Fourier transform has the single 
non-zero coefficient: 

Thus, the first term in the left-hand side of (23) reduces to 

Tr(y^Ay^) - «-iVot,Ayoo, (25) 

where u = U/N^. Similarly, the mean Green function is transitionally invariant in 
space and 'time': Gjj'{T,T') = Gj-j'{T — r'), hence its Fourier transform is diagonal: 

Gl!^'nn'=GlJ^^.'Sr^n', 7 = 0,^. (26) 

Therefore, the trace in the second term of (23) can be written as 

Tr(G^Ay^) = J2 GLAVoo = TriG^AV^o- (27) 

kn 

Substituting (25) and (27) to (23), we come to 

""^^00+ (l + ^)rTrG^ = 0. (28) 

The mean Green function G is related to the mean spin moment (per atom) Sz 
by the formula (see (B.4) in Appendix B) 

Sz = N-^TTrG\ 

Thus, the mean-field equation (28) takes the form 

Vz = -u{l + T])sz, (29) 

Vz = N~^Vqq being the value of the mean exchange field. The conservation of the 
total number of electrons condition {dJ^/ dji = 0) yields the equation on the chemical 
potential (see Appendix B): 

iVe = TTrG. (30) 
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2.5. Equations for the spin fluctuations 

Reduce quadratic form (20) to a sum of squares using the momentum-'frcqucney' 
representation. By Parseval's identity, the first term in (20) can be rewritten as 

^Y^^iAV-r = ^Y1 \^V^mf = u-' E \AV^J'. (31) 

a qma qma 

For the trace in the second term of (20), using (24) and (26) we have 
Tr (GAVGAV) 

kkinni a/37i72 

where 7i,72 = 0,z. Calculations show that the summands with a ^ (3 are equal to 

zero (see Appendix C). Hence 



qma qma 

where q = k — ki, m = n — ni, and 

:)2 p(2) 



l{l+3v)TTv{GAVGAV) = ^ ^V^mX^m^V^^-m = E X^AV^J^ (32) 

qma 

— ni , and 

^'1'" 2 a(Ay«„)5(Ayv„) 

= -^(1 + 3^)tE E GZGlU,n-m Sp{a^^a-a^^a-) (33) 

kn 71 72 

is the unenhanced dynamic susceptibility. Substituting (31) and (32) to (20), we 
obtain 

F(2)(y) = E (u-' - X^m) |AKr„P = E ^qmlAKT^r (34) 

qma ^ ^ qma 

Thus AV^^ are statistically independent Gaussian fluctuations, with the mean squares 
of fluctuations 



3. Local approximation of the SFT 

3.1. Reduction to local fluctuations 

Expressing the mean Green function G in terms of the chemical potential /u, mean 
field Vz and fluctuations (|A]/^p), we come to the closed system of equations (29), 
(30) and (35). However, this system of equations is still excessively difficult for 
practical computations of average quantities, like magnetization. Moreover, as a 
result of the quadratic approximation (34) the fiuctuations at different momenta and 
'frequencies' AV^ become independent, which is an acceptable approximation only 
at low temperatures. Therefore, we proceed with the local approximation of the 
quadratic form (34): 

F(2)(F) = E A^JAV^J-" E^"E l^^cT™!' 

qma a qm 

= E^aA"(i-EiAv,"„r) =E^"^^- 

rv ^ nm. ^ rv 
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Here the coefficient = N^.A"' is related to the mean square of the local fluctuation 
{AV^) by the formula 

\ qm J 

where 

(...)= (/'exp(-^A„Al/Vr)dv) /'...cxp(-^A„AF„VT)dV. 

Taking into account (35), for the mean square of the local fluctuation ('fluctuation', 
for short) we have 

3.2. Summation over momenta and 'frequencies' 

Calculation of the sum (36) follows [9] and yields essentially the same formulae 
but with the renormalization prefactor 1 + 3?? for the susceptibility. By analytic 

continuation from the points icj^ = VlirmT , summation over 'frequencies' is replaced 
by the integration over the energy variable (for details, see [25]): 

where B{e) = {cxp(e/T) — 1)~^ is the Bose function. We discard the temperature- 
independent term with 1/2, assuming that the zero-point fluctuations are already 
taken into account in the initial DOS ^{e) calculated by the density-functional 
method and in the effective interaction constant u. Using the Tailor expansion 
Xq(£) ~ Xq(0) + i<^qe and approximation 

i ^/ T/e, e < eo = (7rV6)T, ^^^^ 



E 



C-^/^-l i 0, £>£o, 

for the Bose function, we come to 

E(|AKrmP) = ;?^-arctan^^, (39) 
qml / 2A«7r 6A« ^ ^ 

m q q 

where 

dlmxq{e) 



A« = 1 - nx"(0), 



de 



(40) 

e=0 



The approximation (38) not only reproduces the behaviour of the Bose function 
-B(e) with respect to thermal energies, but also has the same first moment 
Jq eB{e)d£ = (7rr)^/6, which essentially defines the upper bound eo- Thus, the 
approximation (38) is well justified. Its another advantage is the possibility of 
the straight-forward proceeding to the static limit at high temperatures, when the 
argument of the arctangent in (39) is much larger than unity. 

The function Aq is calculated using the expansion for the static susceptibility: 

Xq(0)«Xo(0) + SV, (41) 
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where Xo(0) static uniform susceptibility, and coefficient B" is obtained from 

the local susceptibility — -^aT^ S^qXq(O)- Substituting (41) in the first equality 

(40), we get 

K = ^0 + (A£ - K)qV¥, (42) 

where Xq = I - uxo{0), = 1 - uxZ{0), and ^ = N'^ Eq9^- The function >p'^ is 
approximated by its mean value: 

V'q^^a-'E^^^L- (43) 

q 

The summation over q is carried out by the integration over the Brillouin zone, 
approximated for simplicity by the sphere of the same volume. Using (39), (42) and 
(43), for the local fluctuation (36) we finally obtain 

where 

al = X^/XZ, bl = il-al)/0.6, Co. = u^l^^T / {QXD . 

Relation (44) is obtained using the integration over the Brillouin zone with the Bose 
distribution and a simple dispersion relation, whose parameters are chosen to give a 
correct value of the local susceptibility. Thus, expression (44) for the moan square of 
the spin fluctuations is self-consistent and does not contain any free parameters. 



3.3. Mean single-site Green function 

To calculate the local susceptibility Xl (-) = XlC*^) +iVL^> replace the mean Green 
function G in (33) by its site-diagonal part: 

^n'n = Gj-j',n ~ GonSjj' = gnSjji , (45) 

so that the Fourier transform is k-indcpondent: Gkn = Sn- Thus, we rewrite (33) as 
xUiuJm) = -^(1 + ^V)TJ2 E fl^ni'^n)5^Mi'^n - icom) Sp(<7TV«<7TV«). 

n 7172 

Using analytic continuation first from points ia;„ and then from iUm, we replace the 
sum over 'frequencies' by the integral over the energy variable: 

x2W = -^(l + 3r?)E / Im(3^^(£)[5^n£-^) 

7172 

+ g^^ie + z)] Sp(a^W^V«))/(£) de, 

where f{e) = [exp((£: — /i)/r) + 1]""'^ is the Fermi function, and 9(e) — g{e — iO). The 
2x2 matrix ^(e) is spin-diagonal: g^{e) = g^{e) = 0, same as Gkn (see (22)). Denoting 
the diagonal elements of 51(e) by g^fie) and gi{s), we have g^{s) = ^[g^{e) + gi{s)] and 
^^(e) = hM^)-9i{£)]- Hencell 

^X(O) ^ _m±M I l^^g^g^) f d,, (46) 
II Recall that Xqm = Xqm due to axial symmetry (for details see Appendix C). 
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=^^^^^^^ flr.g,lrngj-^]ae, (47) 



TT J V ds 

Xm = -^^^^^^ I {lmg',+Imgl)fds, (48) 

= J + (I-.,)^] (-1) de. (49) 

Similarly, in equations for the mean field (29) and chemical potential (30) , we use 
the mean single-site Green function (45) and replace the sum over 'frequencies' by the 
integral over the energy variable. Thus, we rewrite equation (29) as 

V. = -u{l + V)^J Mat - 9i)fde (50) 

and equation (30) as 

= iVa^ J Im{g^ + g^)fde. (51) 

In formulae (46)-(51), the mean single-site Green function is given by 

^-^^^=/ e-.V;-Al(s)-s- ^^' 

where a =t, -i- or ±1 is the spin index, i/(e) is the non-magnetic DOS, and AT,c{e) is 
the fluctuational contribution to the self-energy part, calculated by the formula 

1 + 2aVzg^{e) 

The latter is obtained from the coherent potential equation 

AS = (AT/[1 - ,g(Ay - AS)]-i) 
in the second order with respect to the fluctuations A^ [9]. 

3.4- Correction coefficient 

Final computational formulae of the extended SFT differ from the ones in [9] by the 
renormalization of susceptibility (46)- (49) and the mean field (50) that depend on the 
coefficient rj defined in (21) (with rj — the system of equations reduces to the one 
in [9]). Using the single-site and quasi-static approximations, we come to the following 
expression for ij (see Appendix D): 



2cj {Reg^lmgi+Regilmg^)fde+cj {Reg^lm.g^^+Reg^lmg{) f de 



.(52) 



The impact of the corrections due to the third- and fourth-order terms becomes 
critical at high temperatures. Therefore, coefficient r/ can be estimated by the reduced 
formula 

^7^^/ Re5°Im5Vd£=^C, (53) 

where ^ — (2C^ -|- C^)/3 is the mean fluctuation, and — [g^ + gi)/2. In 
the ferromagnetic region, formula (53) follows from the initial formula (52) in the 
approximation g^ = g^- In the paramagnetic region, formulae (52) and (53) coincide. 
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4. Numerical results 



The extended SFT is investigated by the example of the Invar alloy Feo.65Nio.35. The 
initial non-magnetic DOS 1/(5) (see figure 2) is formed from the two spin-polarized 

densities obtained from the self-consistent calculation for the completely disordered 
Feo.65Nio.35 alloy [26]. The experimental value of the spin magnetic moment per atom 



O 



CO 

(D 
+J 
nS 
+J 

CO 



> 




Figure 2. The DOS of the d band of non-magnetic Fe().(i5Ni().;j5 , obtained 

from [26] { ), and the one smootlied out by convolution with the Loientz 

function of the half-width F = 0.001 { ). The energy e and half-width F are 

in units of the bandwidth W = 9.70 eV. The vertical line indicates the position 
of the Fermi level 6p . 



= 1.70 /iB (where hb is the Bohr magneton), used to determine the effective 
interaction constant u, is taken from [27] . 

Note that we neglect here the fine effects of atomic and/or magnetic short-range 
order (see, e.g., [20, 28] and references therein). Moreover, the magnetic moment 
and DOS u{e) represent the values per averaged atom. However, as stated 
in the introduction, even with those initial data one can calculate the temperature 
dependence of magnetic properties of an alloy in the SFT. 

Figure 3 presents the basic magnetic characteristics for the Feo.65Nio.35 Invar 
calculated within the Gaussian SFT [9]. Clearly, at high temperatures, the calculated 
magnetization curve does not fit well enough the experimental one. For the Curie 
temperature, we obtain Tc = 0.83T^^. But most important, the calculated curve 
m(T) has the inflection (see the discussion in [10]). 

In [13] we took into account the higher terms in the expansion of the free energy 
F{V) by using the simplified expression (53) for the correction coefficient r] with 
c = — O.OISW^^ (W = 9.70 eV is the bandwidth). The calculation gave nearly full 
agreement with experiment for the Curie temperature: Tq = 1.02T^'^'' {T^^ = 520 K 
[27]), for the paramagnetic Curie point: Oc = 1.06T^^, for the effective magnetic 
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Figures. The magnetization m(T)/m(0) ( calculation, odoo experiment [27]), 

the mean square of spin fluctuations (^^ (— ■ ■ —) and (^^ ( ) in units 

of um(0)/(2/iB), the reciprocal paramagnetic susceptibility x~^{T) (— ■ — ) in 

units of r™''//ig, and the local magnetic moment mi^{T)/m{0) ( ) of the 

Invar Feo.csNio.ss calculated in the Gaussian SFT as functions of the reduced 
temperature T/T^"^. 



moment: nies = 0.90m^ ("^eff^ = 3.3 /xb [29]) and for the local magnetic moment 
mi,{T) (see the discussion in [22]). As can be seen from figure 1 in [13], a sharp increase 
of the fluctuations and sharp decrease in magnetization at high temperatures, which 
occurred in [10], disappear in the extended SFT. 

On the whole, the curve for the magnetization in [13] fits the experimental one 
well enough. However, the inflection in the temperature dependence, reported in [10], 
does not vanish entirely. Therefore, in the present paper, we apply the expression (52), 
which alternates with T in a self-consistent way. This finally gives a smooth curve 
without the inflection (figure 4). 

A qualitative analysis of our equations explains the mechanism that leads to a 
jump transition of magnetization in the simple Gaussian theory and its elimination 
in the extended SFT. Indeed, in the presence of the external magnetic fleld h the 
mean-fleld equation (29) takes the form 

V,{h) = -u{l + i^)s,{V,{h) + h). 
Hence for the enhanced magnetic susceptibility we get 

_d£^ ^ Xo ...x 

dh l-u(l + r?)xo' ^ ' 

whore Xo = —dsz/dh is the unenhanced (with constant Vz) susceptibility of non- 
interacting electrons. The fluctuation (44) is proportional to 

1 . , 




The enhanced susceptibility (54) diverges at the critical temperature (the condition of 
the transition from ferro- to paramagnetic state). If, at the same time, the fluctuations 
also increase sharply, i.e. the derivative tends to infinity, then there exist an extra 
(unstable) solution and hence a jump transition to the paramagnetic state. Such a 
scenario takes place in the Gaussian approach (ry = 0), where the susceptibility (54) 
and the amplitude of the fluctuations (55) diverge simultaneously. In the present 
variant of the theory (77 < 0) the fluctuations are weakened, and as a result we observe 
a continuous magnetic transition. 

There are different ways to go beyond the approximation of the non-interacting 
spin fluctuations. In this respect our approach can be treated as a variant of the SCR 
theory. Formulae (54) and (55), together with (53), schematically demonstrate the 
character of our renormalizations. 



5. Conclusions 

We have developed an approximation to the SFT that describes the thermodynamics 
of magnetic characteristics. Our approach takes into account both dynamics and 
non-locality of thermal spin fluctuations, as well as their mode-mode 'frequency' 
interactions. As the initial data, the calculation employs the value of the magnetic 
moment and ab initio DOS, calculated at zero temperature. Further self-consistent 
treatment of thermal fluctuations, including 'large' non-Gaussian fluctuations, makes 
the approximation of the SFT applicable for all temperatures. Particularly, the present 
theory yields a proper second-order phase transition from ferro- to paramagnetic state. 
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Appendix A. Relation between the free energy and Green function 

In this section, we derive the relation between the free energy of electrons in the 
external field and the Green function used in the main text^. Exact expression for 
the constant term of the free energy allows us to obtain the equation on the chemical 
potential as the conservation of number of particles condition in the next Appendix. 

The thermodynamic potential of non-interacting electrons in the external field V 
is defined as 



Oi(y) = -TlnTr 



Trexp(^- H'{V)dT 



(A.1) 



where the Hamiltonian corresponding to the grand canonical ensemble: 

H'iV)=H^ + V{T), 
consists of the Hamiltonian Hq = Hq — i^Nf. and the external field 

Vir) = Vj.„'{r)a%{r)aj„,{T) = 5^Sp(F,(r)p,(r)), (A.2) 

j<y<y' j 
the interaction representation of an operator O is defined as 

Tr stands for the 'time'-ordering operator, and Tr contains the summation over states 
with any number of particles. Note that in H'{V), fi is the chemical potential, 
Ng = aj^(T)aj>(T) is the number of particles operator, and the 
creation/annihilation operators for Wannier states, and pj is the local spin density 
matrix, with elements 

Pjaa' =a^a-aja- (A.3) 

The method to relate the thermodynamic potential ^liV) to the Green function 

G{V) = {z-H\V))-' (A.4) 

is to vary the strength of the external field from to V . To this end, we consider 

H\\)^H'q + \V, (A.5) 

so that -ff'(O) = H'q and -ff'(l) = H'iy). The thermodynamic potential corresponding 
to H'{\) is 

T^exp(^-^ ' H'{X)dT 



= -TlnTr 



Using formula (A.5), we find the derivative of rii(A) with respect to A (cf. [30]): 



F(r)dT) , (A.6) 



dX \ JO IX 

^ The relation between the thermodynamic potential Q,\(y) and Green function G{V), in the case of 
an arbitrary perturbation V , was obtained in [30]. Without a proof the formula, analogous to (A. 12), 
was presented in [31]. 
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where the averaging {. . .)\ of an arbitrary operator O is defined as 
_ Tv((!)r.exp(-^^/^g'(A)dr)) 

General formula (A. 6) holds for any perturbation V, not necessarily one-particle 
operator. For a non-interacting system, substituting (A. 2) for V in (A. 6) and 
rearranging, we write 

l/T 

aOi(A) 



dX 



= T J ^Sp(y,(T)(ft(r))Jdr. (Ai 



The average of the spin density matrix in the interaction representation is related to 
the Green function: 

(p.w(t)), = -(Traj^ir) a+^,{T + 0))^ ^ G^...(t,t). (A.9) 

Substituting (A.9) to (A. 8), we come to 

o \ 

Integration over A between and 1 yields 

f^i(l) -fii(O) = / rTr(VG^)dA, (A.IO) 
Jo 

where Vti{l) = Vti{V) and rJi(O) = -TlnTr cxp(-i7^/r). 

The Green function of non-interacting electrons satisfies the equation 

G^ = Go + \GoVG^, (A.ll) 

where Go corresponds to Hq. Express G^ from equation (A.ll) and substitute to the 
right-hand side of (A.IO). Using the cyclic property of trace, we get 



Oi(l) - f2i(0) = TTr / (1 - AGoy)-^Goy dA 

^0 



= -TTr / — ln(l - XGqV) dA = -TTrln(l - GqV). 
Jo dA 

Using equation (A.ll) and the fact that G^ = G{V), we come to+ 

rJi(y) = -TlnTrexp(-iJ^/T)-TTrlnGo-hTTrlnG(F). (A.12) 

The free energy Fi{V) is related to the thermodynamic potential fii(V') by the 
formula 

rl/T 

Fi{V) = -TlnTr 



T^exp - / H{V)dT 



ni{V) + iiN^ 



where H{V) = Hq + V" is the Hamiltonian corresponding to the canonical ensemble, 
and the number of electrons is fixed. Hence formula (A.12) can be rewritten as 

Fi{V) = -TlnTr exp(-i?o/T) - TTrlnGo +TTrlnG(y). (A. 13) 

Omitting the first and second terms, that do not depend on V, we come to formula 
(4) of the main text. The matrix of the Green function (A.4) for the system with 
electrons reduces to (see formula (5) in the main text) 

GiV) = iz + ^i-Ho-V)-K (A.14) 

+ Recall the formula Trln(j4B) = TrlnA + TrlnB, which is valid for arbitrary matrices A and B. 
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Appendix B. Formulae for the total charge and spin moment 

In this section, we express the mean spin moment and total number of electrons in 
terms of the mean Green function. 

The Green function G{V) is related to the spin density matrix p by formula (A. 9): 

{piT)}v=GiT,T), 

where we write (. . .)v instead of (. . .)i defined by (A. 7). As any Hermitian 2x2 
matrix, the local spin density matrix (A. 3) can be expressed as 

ft 

with the coefficients 

/9^^ = ^SpKp,), ii = 0,x,y,z. (B.l) 

Here a'^ is the 2x2 unity matrix, and a" (a = x, y, z) are the Pauli matrices. Formulae 
(A. 3) and (B.l) yield that the scalar component /9° is equal to one half of the local 
charge: 

1 11 

and the vector component pj is equal to the local spin sj = (sj,s^,s|). Indeed, 
writing the operator Sj in the second-quantisation form, from (A. 3) we find 

Since the spin operator can be represented by the Pauli matrices: s = ^cr, using (B.l) 
we obtain 

s"-^Sp(aV,)=p,". (B.3) 

Relations (A.3), (B.2) and (B.3) lead to 

l{{nj{r))v) = GO(r,r), ((sf (t))^) = (5«(t,t), 

where (. . .) is the averaging over all configurations of the field V with the probability 
density p{V) oc exp(— ri(F)/T). Introduce the Fourier transform 

Gj{t,t') = T^G,-(ia;„,ia;;)e--"-e<-'. 

Using the 'time'-invariance of the mean Green function, we have 
G,-(t, r) = T ^ (iw„) ^ T ^ G,„. 

n n 

Hence, summing over the sites and bands, we come to the expressions for the total 
number of electrons: 

N, = Afd5^((n,» = 2N^tJ2G% = NaTj^SpGjn = TTvG 

j jn jn 

and the total z-projection of spin moment; 

^' = ^dl]((s|)W=^^TrG^ (B.4) 
j 
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Sx and Sy being equal to zero in the ferromagnetic case, since = and = 0. 
Analogously, for the non-interacting electrons we have 

A^e = rTrGo. (B.5) 

Now consider the total free energy T = —TlnZ, where the partition function 
Z is defined by (2). Since jj, is the Lagrange multiplier in the expression for the free 
energy: = Cl + fiN^, the equation on fj, follows from the extremum condition: 

The latter is exactly the conservation of the electrons condition. Since Fq {V) does not 
depend on n, differentiation of with respect to n yields 



0. 



^ - 1 



V 



-Fo{V)/T 



DV{t) 



dFAV) 
djjL 



= 0. 



Using formulae (A. 13) and (A. 14), we come to 

'dF,{vy 



= -rTr(G) +TTrGo = 0. 



Due to (B.5), we finally obtain 

/ dF^jvy 

\ dn , 

The latter is the equation for the chemical potential (30) of the main text. 



-TTrG + N^ = 0. 



Appendix C. Axial symmetry relations 

In this section, we prove the axial symmetry of the magnetic susceptibility in the 
ferromagnetic state: x'^ — (with the 2;-axis chosen along the field), and hence axial 
symmetry of the fluctuations: C,^ = C- 

For given q and m, calculate the susceptibilities 



Aqm — 



1 

2 9(Ay«j5(AFVJ 
2 



kn 7i72 

where 71, 72 = 0, 2. First, using properties of the Pauh matrices, we find 





Xqm 


■yxy 
Aqm 





Xqm — 


Aqm 


XX 

Aqm 













zz 
Aqm 



(C.l) 



where 



V = 
Aqm 



Vqm 



-■^d(l + 3??)r^^Gk„Gk_q_„_m - Gk„Gk_q,„_m)) 
kn 

-A^d(l + 3??)ry^fGk„Gk-q.rt-m + ^kn^k-ci,n-rn) 
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and 

X^L = iA^d(l + 3r,)r^(G°„G^-q,„_„ - Gl„Gi_^^^_Jj . (C.2) 

kn 

Next, we prove that the non-diagonal elements (C.2) of matrix (C.l) are equal to 
zero. Introducing the new indices, we rewrite the second part of the sum (C.2) as 

^kn^k— q, n— m k'+q, — n'+m^— k'— n' ' 

kn k'n' 

Using the property G^^^ = {G"^_^Y of the Fourier transformation (5g„ of the real 
function G'^{t), we have 

^kn^k— q,n— m / ^ V^k^n^ ^k^— q, n^— m/ * 
kn k'n' 

Hence for the non-diagonal element (C.2) we get 



Aqrn 



-27Vd(l + 3r;)T^Im(GLGU 



kn 



Xlm = -^T^E E ^^n<5;:^-q,n-™Sp(a^^a«a^V«), (D.2) 



Since 77 is real, the latter is also real. On the other hand, matrix (C.l) is Hcrmitian, 
hence Xqm must be imaginary. Thus = 0, and the susceptibility x is diagonal in 
the momentum- 'frequency' representation. Formula (36) yields the required relation 
for the fluctuations: ('^ = 

Appendix D. Derivation of formula (52) for the correction coefficient 77 

To calculate the correction coefficient (21): 

7] = (Tr(l))"'^ Tr^GA FGAy ), (D.l) 
where Tr(l) = 2Ni,Nd, we introduce the Gaussian susceptibility 

2 

kn 7172 

i.e. the susceptibility Xqm with no account of high-order terms (jj = 0). Then formula 
(D.l) can be rewritten as 

77 = V Ay" AV" = V y" (lAV" f). (D.3) 

N N T ' I qmAqm'-' ] — q— m N N T ^ A,qm\l"*qml /• V"^'"/ 

^ ^ qma ^ ^ qma 

We approximate (D.2) using the single-site Green function (45) and the quasi-static 
approximation (which implies that the main impact to (D.2) is due to the terms with 
n — m f=s n): 

Xq™ « Xoo = -^^E E ffonffJn Sp [a^V«a^V«] . (D.4) 

n 7172 

Using (D.4), we rewrite (D.3) as 

^ -iVd^E^ooE^l^^cT™!') = "A^E^oor- (D.5) 
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Due to the cyclic property of trace and the anticommutation relations for the Pauli 
matrices: a"a^ = (2(5^^ — l)a^a°', we obtain 

n 

Replacing the summation over 'frequencies' by the integration over the energy variable, 
we come to 

X?o = - v/ + ^^^"^ ~ l)(5^)')/d£. (D.6) 
Substituting (D.6) in (D.5) and using the axial symmetry (^^ = (y), we get 

V = ^ (2CJ Im((5°)2 - {g^') de + cj Mid'f + {ff) de) • 

Recalling that g°{e) = l[gf{e)+gi{e)] and ^^(e) = l[g^{s) - gi{e)], we finally obtain 

^~~f (^^ Im gi + Re 5; Im g^) f de + cj (Re Im + Re g^ Im 5^) / de 

(see formula (52) of the main text) . 
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